† Corresponding author. E-mail:
The flow properties and substrate deposition rate profile, which are the important parameters in electron beam physical vapor deposition, are investigated computationally in this article. Collimators are used to achieve the desired vapor beam and deposition rate profile in some applications. This increases the difficulty measuring boundary conditions and the size of the liquid metal pool inside the collimators. It is accordingly hard to obtain accurate results from numerical calculations. In this article, two-dimensional direct simulation Monte Carlo (DSMC) codes are executed to quantify the influence of uncertainties of boundary conditions and pool sizes. Then, three-dimensional DSMC simulations are established to simulate cerium and neodymium evaporation with the collimator. Experimental and computational results of substrate deposition rate profile are in excellent agreement at various evaporation rates and substrate heights. The results show that the DSMC method can assist in metal evaporation with a collimator.
Electron beam physical vapor deposition (EBPVD) plays an important role in many science and technology areas, such as thermal barrier coatings of turbine blades, the manufacturing of metal matrix composites, supersonic molecular beam technology, large-area corrosion protection coatings, and the enrichment of uranium.[1–3] In EBPVD, a high-energy electron beam is used to heat the metal ingot to create a molten metal pool, where the metal vapor expands from the surface to the vacuum. Collisions between vapor atoms continuously change their trajectories and velocities and then significantly influence their final position and substrate deposition rate, which is of critical importance. Fan et al. researched yttrium vapor deposition in a vacuum chamber both computationally and experimentally.[1] The atomic absorption spectra and deposition thickness profiles obtained by the direct simulation Monte Carlo (DSMC) method and experiments were in excellent agreement at two different evaporation rates. Venkattraman et al. carried out DSMC simulations to research the deposition of copper thin films by adding various molecular models for copper–copper interactions, and a suitable model was determined according to comparisons with published experimental data.[4] Moreover, the influences of thermal nonuniformities at the vapor source were discussed in another article.[5] Balakrishnan et al. investigated the influence of electronic excitation on the flow properties of titanium vapor, and the result compared with experiments showed that the electronic energy is an indispensable factor to consider in natural metal vapor flow.[2] Besides, Thakur and Sahu measured the spatial distribution of several metal vapors experimentally to match Knudsen’s cosine law.[6–8] Chatain et al. measured the velocity, angular distribution, and population density of the ground and first four meta-stable states of gadolinium vapor by time-of-flight, deposition, and laser resonance absorption methods, respectively.[9]
In some applications of metal evaporation (e.g. the enrichment of uranium), a collimator is placed upon the crucible to achieve a desired divergence of the vapor beam and deposition rate profile. As a result, it is difficult to accurately monitor the temperature of the liquid metal pool surface because the collimator shelters it from optical sight.[10] Moreover, the collimator may be heated to make the metal that adheres to the inside walls of the collimator flow back into the crucible to improve the utilization ratio. This means that the sticking coefficient of metal vapor on the inside walls is not exactly 1 in the numerical calculation and experiments. In published papers, scientists usually concentrated on the domain upon the collimator by giving a specific incidence condition at the aperture and ignoring the rest of the domain.[11] However, although good agreements could be obtained at a certain correction, the practical significance of this method is relatively low.
In this article, we carry out DSMC simulations to calculate metal evaporation with a collimator in the whole domain. First, the problem is stated. Then several two-dimensional examples are calculated to quantify the influence of uncertainties of boundary conditions and pool sizes. Next, three-dimensional calculation results are compared with the experimental data of cerium and neodymium. Finally, some conclusions are given.
A schematic diagram of metal evaporation with a collimator is shown in Fig.
First, the surface of the liquid metal pool is divided into an isothermal high-temperature region and a low-temperature region. It is supposed that only the high-temperature region contributes to evaporation according to the formula for the saturated vapor pressure of metal. Then, an approximate area of the high-temperature region, which can be estimated according to the scan pattern of an electron gun, is presented. Next, we give the temperatures and sticking coefficients of the inside walls of the collimator. It is important to note that the values we present here are not necessarily actual. Finally, since the net mass flow rate M can be easily measured experimentally in metal evaporation with a collimator, the temperature of the high-temperature region is adjusted to match the flow rate M at the aperture. Numerical simulation results show that if the net mass flow rate is relatively low, the substrate deposition rate profile would be unrelated to the indeterminate quantities S,
Two-dimensional and three-dimensional DSMC simulations are carried out by employing the variable hard sphere (VHS) model. The collision cross section between two atoms has an inverse power relation with the relative speed in the VHS model. Furthermore, the scattering velocity after collision is isotropic in the center-of-mass system. The relationship between the collision cross section
An orthogonal coordinate system is employed in the simulation. A smaller cell size is adopted near the source surface and aperture to capture the rapid expansion of metal vapor, which is important for accurate calculations. About one million super particles are used. Particles arriving at the substrate are counted to obtain the deposition rate when the evaporation reaches stability.
Two-dimensional DSMC simulations at several net mass flow rates, the results of which are shown in Figs.
The numerical results can be explained as follows. Metal vapor atoms expand to the vacuum collisionally under ordinary circumstances. We suppose a limiting case of weak evaporation, in which the vapor flow above the liquid metal pool reaches the free-molecular or near free-molecular regime. As a result, it can be considered that vapor atoms at the aperture come from only two parts. The vast majority fly from the pool surface directly without any collision while the rest initially collide with the inside wall of the collimator once (the tiny minority of atoms that undergo several collisions is ignored). The profiles of substrate deposition rate shown in Fig.
Figure
A similar structure with a 2 cm×20 cm aperture is used for neodymium evaporation, and the corresponding DSMC simulations are carried out. The deposition weights are measured by weighing the mass changes of patches before and after the experiment at the evaporation rates of 142 mg/min, 183 mg/min, and 266 mg/min, respectively. The deposition rates shown in Fig.
Metal evaporation with a collimator is investigated by the DSMC method. First, some uncertainties of boundary conditions and the size of the liquid metal pool in the simulation are introduced. Then, two-dimensional DSMC calculations are carried out, verifying that the uncertainties have little influence at a low evaporation rate. Finally, three-dimensional DSMC calculations for cerium and neodymium evaporation are performed, and the numerical profiles of deposition rates are compared with those measured by a film thickness gauge and weighting. The results show quantitative agreement under different metals, evaporation rates, and substrate heights. Therefore, it is reasonable to adopt the DSMC method to study metal evaporation with a collimator. It is almost certain that the simulation would be useful for subsequent experiments.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] |